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Abstract 



Rotation-induced time-odd components in the nuclear mean field are an- 
alyzed using the Hartree-Fock cranking approach with effective interactions 
SIII, SkM*, and SkP. Identical dynamical moments are obtained for pairs 
of superdeformed bands 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 151 Tb(l). The cor- 
responding relative alignments strongly depend on which time-odd mean-field 
terms are taken into account in the Hartree-Fock equations. 
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I. INTRODUCTION 



A description of ground and excited states in terms of a mean field is a well estab- 
lished approach in nuclear structure. Most of the nuclear phenomena can be considered 
as manifestations of the mean field properties, or the mean field can be used as a suitable 
framework and the first-order approximation. The mean field results from averaging the 
nucleon-nucleon interactions over states of individual nucleons. The averaging procedure, 
which must take into account the Fermi statistics of nucleons, can be formalized in terms 
of a variational approach, and leads to the well-known Hartree-Fock (HF) selfconsistent 
equations [|TJ. 

Dynamic or time-dependent phenomena can be described by a corresponding time- 
dependent HF method. The nuclear state is then represented by a one-body density matrix, 
which evolves in time according to the Hamilton equations, and represents a motion of a 
wave packet. Such an approach has been applied to genuinely time- dependent problems, 
like nuclear reactions, but it is in fact best suited to describe stationary collective states. 

The nuclear rotation is an example of a collective motion for which linear combinations of 
stationary states of a given spin can be identified with a rotating wave packet. In this case, 
a transition from the time- dependent HF theory to a stationary problem can be achieved 
by introducing rotating intrinsic frame of reference. In this frame, the equations of motion 
are time- independent, however, the resulting density matrix is not invariant with respect to 
the time-reversal operator. As a consequence, the mean field obtained for such a density 
matrix also loses its time invariance, and acquires new terms which are odd with respect to 
the time reversal. 

Properties of nuclear time-even mean fields are known rather well, because they are 
reflected in multiple static phenomena which can be studied experimentally [§]. On the other 
hand, very little is known about properties of the time-odd mean fields. Most studies were 
up to now carried out either within the adiabatic || or semiclassical approximations . In 
particular, the adiabatic approximation to the nuclear translation, rotation, or quadrupole 
motion leads to the well-known Thouless-Valatin || corrections to the mass, moment of 
inertia, or to the vibrational mass parameters, respectively. 

These corrections reflect the fact that a velocity-dependent mean field should be appro- 
priately transformed to the intrinsic frame of reference [0] corresponding to a given collective 
mode. In fact, simple estimations and numerical calculations show that the inertia ob- 
tained for effective forces without the effective- mass terms (m*=m) are very close to those 
for m*<m when the Thouless-Valatin corrections are consistently included. 

Fast nuclear rotation is a phenomenon in which the collective motion should be described 
beyond the adiabatic limit. Such a cranking model has been successfully used in explaining 



numerous high-spin effects in nuclei |10j- In this approach, the properties of the rotating 
mean field explicitly depend on the angular velocity. However, most studies performed up 
to now were done in terms of phenomenological mean fields, which are not selfconsistently 
depending on the rotating states, and, therefore, do not incorporate time-odd terms. Only 
in Ref. |TT[ an attempt has been made to include the effects of a zero-range interaction 



within the Nilsson single-particle mean field. Full selfconsistent cranking calculations [p^-|22 



are still rather scarce, and no explicit analysis of the time-odd mean-field components is 
available. 
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On the other hand, the discovery of the superdeformation [p3|| , and the resulting 



avalanche of very precise data on high-spin states, allow for an attempt to study these un- 
known aspects of the nuclear mean field. In particular, the phenomenon of identical bands 



(see Ref. [24| for a recent review) provides an extremely rich and puzzling information on 
properties of fast rotating states. In terms of the mean field approaches explicitly depending 
on the time-odd components, the calculated identical 7-ray transitions have been obtained 
in two cases only, namely, (i) the HF cranking calculations with Skyrme interaction gave 
yrast band of 194 Hg identical to an excited band in 194 Pb |14| , and (ii) the identical bands in 



152 Dy and 151 Tb (excited band) were obtained in the relativistic mean- field (RMF) theory 
|2l|j . In both studies, the authors invoked the time-odd components of the mean field as a 
crucial element of the obtained results. 

In the present study we aim at a detailed analysis of these time-odd terms in the context 
of the identical bands phenomenon. We focus our attention on two classic and experimentally 
well studied pairs of identical bands, namely, on those in 152 Dy and 151 Tb (excited band), 
and in 151 Tb and 150 Gd (excited band). We perform our calculations in terms of the HF 
cranking method with the Skyrme interactions. 

It was not in our intentions to give here a full account of the time-reversal breaking 
description that the selconsistent HF cranking approach really offers and the experiment 
may test. Such a field of studies seems to emerge now in relation with the fast progress in 
the available experimental information. In particular, our choice of the illustrative material 
purposely minimizes rather than maximizes a possible magnitude of the time-odd effects. 
This is so, because we mainly discuss the role of the 7r [301] 1/2 Nilsson orbital known to 
interact only weakly via the Coriolis term in the hamiltonian. 

In all standard approaches, such as the Woods-Saxon (WS), Nilsson, HF, or RMF crank- 
ing models, one signature member of this orbital is to a good approximation given by a 
staight line as a function of the angular velocity. This type of dependence can be opposed 
to dramatic changes of other aligning or interacting orbitals which are much more sensitive 
to the rotating field and hence to its time-odd components. On the other hand, the specific 
behavior of the discussed orbital induces a regular behavior of various observables, and has 
a great advantage in the fact that the level repulsion or level crosings do not disturb the 
comparisons of primary interest here. 

In Section y we discuss the Skyrme interaction and energy density with a particular 
emphasis on the relations between the time-even and time-odd components of the mean 
field. The HF cranking calculations are presented in Section |TJ, where several variants of 
the Skyrme functional are used and the influence of the time-odd terms on the rotational 
properties is analyzed, and conclusions are presented in Section [TV]. 



II. SKYRME ENERGY DENSITY 



The Skyrme energy functional [[25|,[2(| is a three-dimensional integral 

S = J d 3 rH(r) 



(2.1) 



of the energy density 7i(r) which can be represented as a sum of the kinetic energy and of 
the potential-energy isoscalar and isovector terms 
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where 



H(r) = ^-T + H (r)+H 1 (r), (2.2) 



Utir) = nr n (r) + nf d (r), (2.3) 



with 



nr n (r) = c t p P 2 t + cpfHApt + c; Pt r t + c t J j 2 t + c? J Pt v-j t , 

Hf d (r) = Qs 2 t + C t As s t -As t + Cjst-Tt + C J t j 2 + c7 J s t -(Vxj 



(2.4) 



th 



and the isospin index t can have values or 1. Coupling constants C t have superscripts 
corresponding to various density functions appearing in Eq. (|2.4|). All of these coupling 



constants can, in principle, depend on particle densities, see Refs. pT| , |2"8|j , but the most 
common choice |^5[^9j restricts the density dependence to the Cf and Cf terms. In Ap- 
pendix [A] we give coupling constants C t expressed through the traditional parameters of the 
Skyrme interaction. 

Apart from the density-dependent coupling constants, the isoscalar energy density, 
Ho(r), depends on the isoscalar density functions and the isovector one, H\(r), depends 
on the isovector density functions. For the particle densities, p t , the isoscalar and isovector 
parts are defined in the usual way as a sum and a difference of the proton and neutron 
contributions, respectively, 

Po = Pn + P P , Pi = Pn - P P , (2.5) 

and analogous expressions are used to define other densities. 

Altogether we have to consider six position-dependent density functions p, r, j, s, T 



and J with definitions given in Ref. p6[ . Whenever the isospin indices are omitted, we 
understand that the symbols may refer either to the isoscalar (t=0) or isovector (t=l) parts. 

In the energy densities ( |2.4|) there appear two scalar time-even density functions, p 
and r, one vector time-odd, j, and two pseudo- vector time-odd ones, s and T, and one 
pseudo-tensor time-even density function, J. The vector time-even density J which ap- 
pears in Eq. ( |2.4j ) is given by the antisymmetric part of the pseudo-tensor density, i.e., 
J\=Huv £\uvJuv, and is not an independent quantity. Terms in Hf cn (r) and 7-^° dd (r) are 
bilinear in time-even and time-odd densities, respectively. Therefore, we denote them by the 
superscripts "even" and "odd". This notation is only used to indicate the dependence on 
two different classes of densities. It should not be confused with the fact that both 7^ ven (r) 
and 7i° Ad {r) are of course even with respect to the time reversal. To every term in 7ij ven (r) 



there corresponds an analogous term in TC° (r), as seen in Eq. 

A. Mean fields 

By varying the energy density (|2.4| ) with respect to the six density functions p, r, j, s, 



T and J one obtains the mean fields. Details of calculations are presented in Ref. |26]. Here 
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we only repeat the final results conforming to the notation introduced above. Time-even 
and time-odd mean fields are obtained by a variation of 7if en (r) and 7i° dd (r), respectively, 
and they read 



(2.6) 



rr n = _ v . M ,(r)V + U t {r) + ± [V a • B t (r) + B t (r) ■ Va 

Tf d = -V • (a ■ C f (r))V + o-.S t (r) + ±(V • I t (r) + I t (r) ■ V). 



These fields are given by the six potential functions U, M, I, E, C and B which have tensor 
transformation properties respectively identical to the six density functions on which they 
depend through the following formulae 

U t = 2C p tPt + 2C t Ap A Pt + C T t r t + C t VJ V ■ J t , (2.7a) 

S t = 2C s t s t + 2C? s As t + CjT t + C? j V x j t , (2.7b) 

M t = C T tPu (2.7c) 

C t = Cfs t , (2.7d) 

% = 2C(j t - C? J Vp t , (2.7e) 

J 4 = 2C{jt + C t Vi V x (2.7f) 



The tensor gradient operators in Eqs. (|2.6|) and (|2.7e|) are defined [^] as (Ver) MI/ =V M er,/ and 

Due to the density dependence of the coupling constants Cf and C t s one has to add to 
the isoscalar potential energy Uq the rearrangement terms [JTJ resulting from the variation of 
these coupling constants with respect to the isoscalar particle density, i.e., 

% = H( ^ + ^4] ■ (2-8) 



d Po rt dp. 



o 



These terms introduce an explicit dependence of the complete time-even isoscalar potential 
Uq+Uq on the isovector density pi and on the time-odd densities sj. Similar terms have to 
be also consistently taken into account whenever some other coupling constants are assumed 
to be density dependent. 

Finally, the neutron and proton hamiltonians, h n and h p , are obtained by combining the 
kinetic energy with the isoscalar and isovector mean fields: 



/, h 2 a I -pevcn I "podd i pcvcn i podd 

' ln ~ 2ni "r i 0"r i 0+ i l+ i l' ^9) 

z, h 2 A I peven _i_ podd peven podd 

n P ~ 2m a " hi +1 1 1 1 1 ■ 

We are now in a position to discuss different time-odd terms in the mean fields. It is clear 
that the time-odd mean fields, T° dd , in Eq. (|2.6|) directly result from the "odd" part, 7i° dd (r), 
of the energy density in (|2.4j), and depend on 10 time-odd coupling constants Cf, C t As , Cf, 
C{, and for t=0 and t—1. Similarly, the time-even mean fields depend on 10 time-even 
coupling constants Cf, C Ap , CJ, C/, and Cf J . For the Skyrme interaction, the time-odd 
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coupling constants are linear combinations of the time-even ones (see Appendix |AT), and 
therefore the time-odd mean fields are uniquely determined from the time-even mean fields. 
Since the time-even fields are tested against numerous experimental observations which have 
a static character, they are much better known then the time-odd ones. In fact, the Skyrme 
force parameters have been in the past almost uniquely fitted to the static properties only. A 
description of dynamic properties which do depend on the time-odd fields does not therefore 
(for the Skyrme force) require new parameters to be introduced and fitted. 

On the other hand, one sometimes adopts a different point of view by considering the 
energy density to be a more fundamental construction than the Skyrme interaction itself 
f30|j . In such a case, all 20 coupling constants of Eq. (|2.4|) should be treated and adjusted 
independently. However, in the next Section we show that some relations between time-odd 
and time-even coupling constants have an origin in the local gauge invariance of the energy 
density. 



B. Local gauge invariance 



As noted in Ref. |[26|| , in the energy density derived from the Skyrme interaction the 
kinetic density r and the current density j appear in the characteristic combination of 
Pt T t~3t- The same is true for two other pairs of densities appearing together in the com- 
binations SrT t —j\ and (p t 'V-Jt+St ■ (Vxj t )). This gives the following relations between 
three pairs of time-even and time-odd coupling constants 

C{ = -C[, (2.10a) 
C/ = -Cf, (2.10b) 
C? j = +C? J . (2.10c) 

The Skyrme functional has now the form 

H t (r) = C p tP \ + C s t s 2 t + C^p t A P t + C? s s t ■ As t 
+ Cl( Pt T t -fi)+C?(s t -T t -J*) 

+ C7 J ( Pt V -Jt + SfiVx j t )). (2.11) 

In Ref. [^fj this structure has been interpreted as a result of the Galilean invariance of 
the Skyrme interaction. However, it has a deeper origin in the fact that the Skyrme force is 
locally gauge invariant. In order to illustrate the role of the locally gauge invariant, velocity 
dependent interactions let us consider an arbitrary finite-range and non-local, but velocity 
independent interaction given by 

V = V{r' i y 2 -r ll r 2 ). (2.12) 

To simplify the notation we disregard for a moment the spin and isospin variables. V 
describes an interaction process where the particles 1 and 2 are located at r\ and r 2 before 
the interaction, and at r[ and r' 2 after the interaction. When the system of particles is 
described by a one-body density matrix p(r, r'), its HF interaction energy reads 
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£ int = J dr' 1 dr' 2 dr l dr 2 V{r' 11 r' 2 ]r l ,r 2 ) x 

(p(r 1; r[)p(r 2 , r' 2 ) - p(r 2 , r'Xn, r' 2 )). (2.13) 
For a local gauge transformation of the many-body HF wave function |\&), 

|tt')=exp^y>(r.,-U|*>, (2.14) 



3=1 

where 4>(r) is an arbitrary real function of the position r, one obtains the following gauge- 
transformed one-body density matrix 



p'(r, r') = exp 



{i{<j>{r)-<t>{r'))}p{r,r'). (2.15) 



The general interaction energy ( 2.13 ) is not invariant with respect to such a transforma- 
tion. However, when the interaction is local 0], 

V(r[, r' 2 ; n, r 2 ) = S^-rMr^-r^V^, r 2 ), (2.16) 
the corresponding interaction energy, 

£ int = J dr 1 dr 2 V(r 1 ,r 2 ) x 

[p{n)p{r 2 ) - p{r 2 ,n)p{ri,r 2 )), (2.17) 



becomes invariant with respect to the local gauge. The direct term is invariant because 
it depends only on the gauge-invariant local densities, which are denoted by the single 
argument, i.e., p(r)=p(r,r). On the other hand, in the exchange term the gauge factors 
coming from two density matrices cancel one another. 

When the density functions defining the Skyrme energy functional ( |2.4j) are calculated 
for the gauge-transformed density matrix ( j2.15| ) one obtains the following relations: 

p' = p (2.18a) 

T > = T + 2j ■ V0 + p(V0) 2 , (2.18b) 

4 = s k , (2.18c) 

fk=h + pVk<f>, (2.18d) 

^ = T fc + 2^J M V ; 0,+ Sfc (V0) 2 , (2.18e) 
i 

J' M = J H + Si V fc 0, (2.18f) 

and the three characteristic combinations of density functions, which appear in the energy 
density of Eq. (|2.11| ), are then explicitly gauge invariant. Transformation properties of r 
and j allow to interpret ClS db velocity field, 

v = (2.19) 

m 

which shows that the flow of matter obtained through the gauge transformation is irrota- 
tional, Vxt7=0. 

The local gauge invariance of the Skyrme interaction reflects the fact that its velocity 
dependence has been introduced only to simulate the finite range effects of the effective 
interaction. In this way the Skyrme interaction conserves the local gauge invariance of a 
velocity-independent finite-range local interaction, such as the Gogny force M , for example. 
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1. Translational motion 



The Galilean invariance is a special case of the local gauge invariance, for which the 
phase in Eq. (|2.14 ) is given by 



<j>(r) = ^f, (2.20) 

where p is a constant linear momentum of the boost transformation. Since the interaction 
energy does not depend on p, the total energy (|2.1|) is only modified through the kinetic 
energy [the first term in Eq. (2^2)]. Using the transformation property of r, Eq. ( 2.18 b| ) , we 



have the energy increase under the boost transformation, 

A £b00St = / 2 

2m 

equal to the translational energy of the boosted system. This result holds for an initially 
stationary solution (i.e., for vanishing currents, j=0), however, due to the transformation 
property of j, Eq. ( |2.18d| ), the boost transformations can be added to one another by 



adding the corresponding momenta p. For the boost transformation, one obtains the obvious 
velocity field (|2~T91 ), i.e., 

^boost = P ^.22) 
m 



2. Rotational motion 

Since the velocity field ( |2.19| ) obtained through a gauge transformation is irrotational, 
it cannot correctly describe physical rotations of nuclei. This is so, because the nuclei 
basically rotate as rigid bodies (at least in the independent-particle approximation ||), and 
the velocity field of a rigid-body rotation, v Tlgld =u} x r, has a non-zero curl, Vxi? rigld =2ct;, 
i.e., is not irrotational. Of course, interactions and nucleon-nucleon correlations (pairing) 
introduce an irrotational component in the velocity field (moment of inertia decreases below 
the rigid-body value), but this field is never entirely irrotational. 

In an analogy to the boost transformation, one may try to induce the rotation of a 
many-fermion system by adding for all particles a constant value, j x , to their projections of 
the angular momentum on a fixed (say x) axis. Such a procedure can be realized in terms 
of the twirl transformation given by the gauge function 

m = j,arctan(^ ^ 

Its velocity field has the form 

v twiTl = ^ — -u> x r, (2.24) 

murj(y, z) 1 

where uj is the vector of angular velocity oriented along the x axis, and r](y, z) is the distance 



to this axis. One can see that the velocity field (|2.24| ) is singular at the rotation axis. The 



S 



energy increase resulting from the kinetic energy density (|2.18b ) is then infinite for any 
arbitrarily small j x . The irrotational velocity field ( 2.24 ) is of course very much different 
from the rigid-rotation velocity field, i> ngld , even though both contain the same factor uj x r. 

This illustrative example shows that the nuclear rotation cannot be realized by demo- 
cratically distributing the angular momentum among all particles. In the cranking approx- 
imation, at a given value of the angular velocity some particles receive larger contributions 
(aligning states) and some smaller contributions (high-.fr states). A precise distribution 
cannot be found without actually solving the quantal cranking equations. This example 
also shows that the gauge-invariant interaction must contribute to the rotational energy, 
contrary to what happens in the case of the translational motion. 



III. HARTREE-FOCK CRANKING CALCULATIONS 

In the present study we have performed the Hartree-Fock (HF) calculations of super- 
deformed rotational bands using the Skyrme effective interactions. The calculations have 
been done using the numerical code HFODD which employs a three-dimensional cartesian 
deformed harmonic oscillator (HO) basis to describe the single-particle wave functions. The 
details concerning the HFODD code will be presented elsewhere here we only give a 
few of its basic parameters pertaining to the present application. 

The calculations have been performed using a fixed basis given by the HO frequen- 
cies 7iu;_l=11.200 and foa;j|=6.246 MeV in the directions perpendicular and parallel to the 
symmetry axis, respectively. These values have been obtained by standard prescriptions 
developed for diagonalizing the deformed WS hamiltonian in the HO basis, and cor- 

respond to the WS potential with deformations /3 2 =0.61 and /3 4 =0.10. The basis has been 
restricted to a fixed number, M, of basis states having the lowest single-particle HO ener- 
gies €Ho = ( n x+n y +l)hu)j_+(n z +7})hu>\\. The actual calculations have been performed with 
M=306. This corresponds to the maximal numbers of oscillator quanta equal to 8 and 15 
in the perpendicular and parallel directions, respectively. 

The stability of results with respect to increasing the size of the HO basis has been 
tested by performing calculations with M=604, which introduces basis states up to 11 and 
20 quanta in these two directions. It has been found that the rotational characteristics 
of the studied nuclei are almost independent of such an increase. For example, numerical 
inaccuracies in the dynamical moment and in the total angular momentum I can be 
estimated to be smaller than 0.2 fr 2 MeV _1 and 0.1 ft, respectively. Inaccuracies of relative 
values between different angular frequencies lu, or between different nuclei, are smaller than 
these estimates, because the numerical errors then cancel out. 

In the present study we aim at investigating the role of different time-odd terms in the 
selfconsistent mean fields Q2.6|) obtained for rotating superdeformed nuclei. As discussed in 
Section [T|, this can be done by considering different values of 10 time-odd coupling constants 
appearing in the Skyrme energy density, Eq. fl2.4|) . For every given set of values of coupling 
constants we perform full selfconsistent calculations within the HF cranking method. In the 
present study, pairing correlations are not taken into account. Below we separately discuss 
three cases corresponding to (i) the complete Skyrme functional, (ii) the Skyrme functional 
with certain time-odd terms omitted, but with the gauge invariance preserved, and (iii) with 
the gauge invariance violated. 
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As discussed in Section |I B| > the energy functionals corresponding to the standard Skyrme 
interactions preserve the local gauge invariance. The 20 coupling constants appearing in 
Eq. ( p.4|) are then restricted by 6 conditions ( |2.1Qa| )-( |2.1Qc|) , which leads to the energy 
density of Eq. (|2.11|) . In Table |, the remaining 14 coupling constants are listed for the SIII 
|29| , SkM* f54| , and SkP |35| Skyrme interactions. The values of the density-dependent 
coupling constants Cf and Cf are given for the vacuum (p=0) and for the nuclear-matter 
saturation density (p=p nm ) characterizing a given force. 

The time-even coupling constants corresponding to the SIII, SkM*, and SkP interac- 
tions are rather similar. The main difference consists in different values of the isoscalar- 
effective-mass coupling constant Cq which is equal to for SkP (effective mass m*/m=l), 
and 44.4MeV-fm 5 and 34.7MeV-fm 5 for SIII and SkM*, respectively (effective masses 
m*/m=0.76 and 0.79). Apart from that, the absolute values of the Cf coupling constants are 



larger for SkP than for SIII and SkM*, which gives better symmetry-energy properties |36 



within the SkP parametrization as compared to the other two forces. On the other hand, the 
SkM* parameters have been adjusted so as to properly describe the surface energy at large 
deformation, and therefore this interaction was successfully applied in numerous studies of 
sup er deformation . 

Values of the time-odd coupling constants corresponding to the three Skyrme interac- 
tions, Table |, differ much more than those of the time-even ones. This illustrates uncer- 
tainties in determining the time-odd components of the mean field. The differences partly 
result from the fact that SkP has been adjusted to give attractive matrix elements in the 
pairing channel, while those given by SIII and SkM* are repulsive. One should bear in 
mind, however, that for the Skyrme interaction the time-odd coupling constants are unique 
functions of the time-even ones (see Appendix |S|) and, in principle, one has no freedom for 
an independent readjustment. Such a readjustment is possible only if we consider the HF 
theory based directly on the energy density functional and not on the Skyrme interaction. 

In the previous applications of the Skyrme force to nuclear rotation JT2|--|il)|] , in the 
energy density (|2.11|) the terms C t As and Cf have been neglected in order to facilitate the 
calculations. The first of these terms gives purely time-odd contribution to the mean field, 
while the second one gives both time-even and time-odd contributions, because the gauge 
invariance implies that C 4 T =C/, Eq. ( |2 . 1 0b|) . In fact, the term C( has also been usually 
neglected in most Skyrme parametrizations applied to problems where the time-reversal 
symmetry is concerved, cf . Ref . |H| . Since the HFODD code is organized in a different way, 
omitting some terms would not provide any serious simplifications, and the code may in fact 
handle the complete Skyrme functional. This gives us a possibility to test the importance 
of different terms for the rotational properties of nuclei. 

There exist several high-spin observables which, when calculated from the HF solutions 
(wave functions) behave differently depending on whether various time-odd terms are in- 
cluded or not. This offers, in principle, a possibility of both better readjustment of the in- 
teraction coupling constants and better understanding of the underlying mechanisms. Some 
of the physical quantities, as e.g. alignments and j7^ 2 ^-moments of intruder orbitals, are well 
recognized as responding strongly to the Coriolis and centrifugal interaction effects. The 
same quantities are expected to respond relatively strongly to the time-odd terms in the 
mean-field hamiltonian. Similarly, various families of orbitals having large high-j compo- 
nents are systematically responsible for such precisely measurable effects as band-crossings 
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(back- or up-bending effects) and related phenomena. 

In this study we would like to illustrate the effects of the time-odd terms on yet another 
seemingly more subtle mechanism related to identical bands, leaving the aformentioned 
analysis of intruder orbitals for a later investigation. 

In the following Sections we present results of the HF cranking calculations for the 
yrast superdeformed bands in dysprosium, 152 Dy(l), and terbium, 151 Tb(l), and for the 
first excited bands in the corresponding isotones, 151 Tb(2) and 150 Gd(2). According to the 
standard notation, the numbers in parentheses refer to numbers attributed in experimental 
studies in connection with relative intensities of gamma transitions. The experimental data 
are taken from Refs. |3?H39f where the most recent and precise results are given. 

The pairs of superdeformed bands, 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 151 Tb(l), are iden- 
tical with a very high precision pOfl . For each pair, the 7-ray transition energies E 1 are 
identical up to 2keV. The identity of the bands can be characterized in two ways, (i) by 
their relative alignments, and (ii) by their relative dynamical moments. The relative align- 
ment SI is defined as a difference of spins in two bands at fixed angular velocity u. Similarly, 
the relative dynamical moment is the difference of at the same value of uj. In 

calculations, the latter is a derivative of the former with respect to the angular velocity. 

In the present study we have fixed the yrast proton configurations of 152 Dy and 151 Tb to 
(16,16,17,17) and (15,16,17,17), respectively. The numbers in parentheses denote the num- 
bers of lowest states occupied in the parity-signature blocks ordered as (++,H — , — h, ), 

where the signs denote the parity quantum number and the sign of the (imaginary) signature 
II quantum number, r=±i. One somtimes uses a notation based on the signature index a 
which equals to — l/2(+l/2) for r=+i(—i). For all bands considered in the present study, 
the neutron configurations are fixed at (22,22,21,21). The yrast configuration in 151 Tb corre- 
sponds to a hole in the 16-th orbital in the ++ block, 16 ++ , which in the standard notation 
is described as the 7r[651]3/2 (r—+i) or 1164 Nilsson intruder orbital. 

Since in the experiment the angular velocity is associated with half of the transition 
energy E 7 , and the above pairs of bands have identical transition energies, the relative 
alignments must be close to a half-integer value. This is simply a consequence of the fact that 
the spins in the odd and even nuclei are half-integer and integer, respectively. Departures 
from half-integer relative alignments can only be caused by differences in the 7-ray transition 
energies, which are very small for the two pairs of bands studied here. Since the values of 
spins have not yet been measured, the relative alignments are known up to an additive 
integer value and a theoretical input is necessary if one wishes to put forward one value or 
another. 

Already at a very early stage of the identical band studies, it has been suggested HI 



that the first excited bands 151 Tb(2) and 150 Gd(2) correspond to the 7r[301]l/2 (r—+i) holes 
(signature index a=— 1/2) in the yrast states of the respective cores (see also Ref. [fffil ). In 
the present study we have followed this interpretation and we have constructed the excited 
bands by creating a hole in the 17-th orbital of the — h block, 17 h . 

In all the WS and HF cranking calculations, the single-particle routhian 7r [301] 1/2 (r—+i) 
increases with rotational frequency with a constant slope of about +0.5 fr, i.e., it has the 
single-particle alignment (the average value of the projection of angular momentum on the 
rotation axis) close to — 0.5 fo. Therefore, a hole created in this orbital must lead to relative 
alignment of about +0.5/1. In this paper we adopt this half-integer value for fixing the 
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unknown integer additive constant, 5Io=+0.5h, required to extract relative alignment from 
experimental data. Within this choice, the experimental average relative alignment for the 
151 Tb(2)- 152 Dy(l) pair of bands equals to +0.564(18) H. The error given here is the average 
error resulting from experimental errors of transition energies. The experimental relative 
alignment for the 150 Gd(2)- 151 Tb(l) pair increases slightly with the angular velocity and 
has the average value of +0.479(14) h. In fact, experimentally, only the average departures 
from the half-integer constant 5I , (5I)=6I +0.064(18)h and (<H}=Wo-0.021(14)7i, are 
established. 

A. Complete Skyrme functionals 

In Figs. and @ we show results of calculations for the 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 
151 Tb(l) pairs of bands, respectively. The complete Skyrme energy-density functionals of 
the SIII, SkM*, and SkP forces have been used. Parts (a) of these Figures (bottom) present 
the dynamical moments as functions of the angular velocity uj. 

1. Dynamical moments 

For 152 Dy(l), Fig. [I], all forces overestimate the experimental values of by about 
5-10%. All three forces give very similar results, within 2fr 2 MeV~ 1 . A similarity of results 
obtained for different forces is also visible in Table [TT], where we give the values of proton 
quadrupole moments calculated at o;=0.5 MeV/h. The SIII interaction gives values larger by 
only 0.3-0.4 b as compared to those given by SkM*, while the SkP force leads to intermediate 
results. A similarity of results for changes of proton quadrupole moments induced by creating 
holes in proton orbitals is even more pronounced. A polarization by the 7r[301]l/2 (r—+i) 
hole gives, for different forces, the changes between 0.12 b and 0.15 b. That induced by the 
7r[651]3/2 (r—+i) hole gives values between —0.94b and —1.05b. The HF proton quadrupole 
moment of 152 Dy(l) agrees very well with the result obtained within the Nilsson single- 
particle potential p2| . 

In Fig. 0(a) we show the dynamical moments calculated for the 151 Tb(l) band. 
Due to the hole in the 7r[651]3/2 (r=+i) intruder orbital, decreases here with to much 
faster than that of the 152 Dy(l) band. However, at c<j~0.45 MeV/h one obtains values of 
j( 2 )~94^ 2 MeV- 1 which are almost identical for all the forces and for both bands. This 
contradicts simple perturbative estimates. Indeed, a hole in the 152 Dy(l) core should in 
principle cause a decrease of the moment of inertia due to the smaller mass and deformation 
of 151 Tb(l), and also due to the fact that the 7r[651]3/2 (r—+i) routhian has negative second 
derivative with respect to u and, therefore, the hole in this orbital should bring negative 
perturbative contribution to . Nevertheless, at cj~0.45 MeV/h the values calculated 
for 152 Dy(l) and 151 Tb(l) are the same. This illustrates the fact that the polarization effects 
obtained by selfconsistent calculations do not necessarily follow perturbative estimates. 

In 151 Tb(l), the SIII and SkM* forces give a better agreement with data than SkP. This 
can be attributed to a different effective masses, m*/m=l for SkP and 0.76-0.79 for SIII and 
SkM*, which leads to different time-odd components in the mean field (see the next Section). 
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Then, the interaction between the intruder orbital 7r[651]3/2 (r=+i) and the time-odd mean 
field is modified, and gives a more significant departure from experiment. 



2. Relative dynamical moments 

In parts (b) of Figs. [1] and we present the relative dynamical moments SJ 1 ^ calculated 
for the pairs of bands 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 151 Tb(l), respectively. One should 
note that the scale in (b) is enlarged five times with respect to that in (a). At u;>0.3 MeV/h, 
both pairs of bands have dynamical moments identical up to about 1 ft 2 MeV _1 , with lighter 
isotones having slightly larger values. In the scale of (a) this would lead to curves identical 
up to the size of the data marks. 



This result confirms the observation [14] that the identity of dynamical moments can 



in fact be obtained in self-consistent theories, and does not necessarily follow semiclassical 
|43| or perturbative estimates. The experimental values of SJ^ are not shown in parts 
(b), because they would be scattered between ±1.5 fr 2 MeV _1 with errors between 1 and 
2.5 fc 2 MeV _1 , i.e., in the scale of (b) they would cover the whole presented region of SJ^ 2 \ 

3. Relative alignments 51 

The calculated relative alignments SI, shown in parts (c) of Figs. |T| and El do not re- 
produce experimental results with a sufficient precision. For the 151 Tb(2)- 152 Dy(l) pair 
[Fig. 0(c)], in the region of u where the data are available, one obtains a gradual increase 
of SI by almost 0.5 h. For the SIII and SkP forces the value obtained at u=0.3 MeV/h is 
correct, but for SkM* the whole curve is additionally shifted up by about 0.5/1. For the 
150 Gd(2)- 151 Tb(l) calculations this increase becomes smaller (SIII and SkP) and similar to 
the small increase seen in the experimental data. However, the values of SI are still slightly 
(SIII and SkP) or significantly (SkM*) too large. 

Values of relative alignments can be translated into the differences SE y of 7-transition 
energies between the two bands. Using a linear local expansion of spin as function of the 
angular frequency we obtain that 



r j £ j 

SE 1 ~ 2%Su ~ 2h ^ {2) ° , (3.1) 



where Suo is the difference of frequencies at spins of physical states in two nuclei. Hence 
the departures of calculated relative alignments SI from SIo=0.5h correspond to the values 
of 5Ey between and 10 keV, while the measured values are between and 2 keV. In the 
following two sections we present calculations obtained with modified Skyrme functionals. 
In this way we try to analyze the influence of the time-odd terms on the relative alignments 
discussed here. 



B. Modified gauge- invariant Skyrme functionals 



As discussed above, the gauge-invariance conditions ( |2.10a| )- (|2.10c|) restrict values of 6 



time-odd coupling constants and leave a freedom to modify the values of C? and C 4 As . If one 
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decides to depart from the complete Skyrme-interaction functional, in which the time-odd 
coupling constants are uniquely defined by the time-even ones (Appendix [A]), one may, in 
principle, use arbitrary values of Cf and C As . However, independent variations of these 
coupling constants have never been considered in the literature, and their effects are, up to 
now, unknown. Therefore, in the present study we restrict our analysis to the functionals 
in which one or more of the SkM* coupling constants, Table H, are assumed to be equal 
to zero. Moreover, in order to further restrict the number of possible variants, we only 
consider simultaneous modifications of the isoscalar and isovector coupling constants of a 
given species. In the frame of gauge-invariant functionals, this leaves us a possibility of 
putting Cf=0 and/or C As =0. 

Another modified gauge-invariant functional can be obtained by removing the term 
s t -T t —J 2 , i.e., by putting C t T =C/=0, in accordance with Eq. ( p.lObJ ). This leads to a modi- 
fication of time-odd and time-even mean fields. However, the term Jf was anyhow neglected 
in most parametrizations o the Skyrme forces used for time-even studies, and in particular 
in SkM*. Therefore, below we discuss four possibilities corresponding to (i) the complete 
functional, (ii) C t T =C/=0, (iii) C A *=Cf=C/=0, and then (iv) Q=C As =C t T =C/=0. The 
results are presented in Figs. ^| and ^ for the pairs of bands 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 
151 Tb(l), respectively. We have studied several other possibilities, but principal conclusions 
can be drawn from these four cases. 

For 152 Dy(l), the dynamical moments , Fig. |3|(a), are (at high spin) sensitive only 
to the C t T =C/ coupling constants. The influence of C t As and CI is visible only at low spins, 
and moreover, their effects have opposite signs and partially cancel one another. Removing 
from the functional the term s t -T t —J 2 decreases the dynamical moments, and, at the high- 
spin-end of the band, brings the calculated value down to the experimental result. At lower 
spins one obtains a smaller than previously overestimation of the data (by about 5%). 

A very interesting results is obtained for the relative dynamical moments 5J^ 2 \ Fig. 0(b). 
For any studied combination of the time-odd mean fields taken into account, the dynamical 
moments in both nuclei are almost identical. Even if the values of J 1 ^ depend on whether 
the term s t -T t —J 2 is taken into account or not, the differences of J^ 2 ' do not depend on it 
at all. This is a rather general observation, valid also for other cases of modified functionals 
discussed below. It means that the identity of dynamical moments is fairly independent of 
at least some of the details of effective interactions, and can probably be attributed to rather 
simple geometric effects. 

The relative alignments 51, Fig. |||(c), very strongly depend on the presence of the 
s t -T t —J 2 term and are almost independent of the s 2 and s t -As t terms. As soon as 
C i T =C/=0, the value of 51 at u;=0.3 MeV/h comes down to the experimental result. On 
the other hand, the increase of 51 as function of the angular frequency is still too large, and 
does not give the experimental identity of the 7-transition energies, Eq. (|3.1| ). 

In Fig. [|, for the pairs of bands 150 Gd(2)- 151 Tb(l) we present a similar analysis of the 
role of different time-odd terms in the Skyrme functional. In this case, the influence of 
the s t -T t -J 2 t term on J® is weaker than for the 151 Tb(2)- 152 Dy(l) pair, while its removal 
leads to a good agreement with experiment, with only a slightly too slow decrease of 
as function of the angular velocity. At the same time the relative alignments become much 
closer to experimental data, but an overestimation by a few tenth of H persists at all angular 
frequencies. The influence of the other two time-odd terms considered here, s 2 and s t -As t , 
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is much weaker, although a non-negligible influence on relative alignments can be seen. 

We conclude this section by noting that the importance of the gauge-invariant time-odd 
terms for the rotational properties is not simply related to the order of the given term. In 
fact, the zero-order term , which depends on the density of spin, has a very small influence 
on the results. At the same time, two second-order terms, s t -As t and s t -T t —J^, which in the 
previous studies have been simultaneously neglected, have a small and rather large influence, 
respectively. On the other hand, the latter term, which in principle cannot be neglected 
because of its magnitude, leads to larger deviations from experimental data, as compared 
with calculations which do disregard it. Of course, the present analysis is not sufficient for a 
more accurate derivation of the magnitude of coupling constants from experiment. This can 
only be done by a simultaneous consideration of many different available data, and moreover, 
should also involve a careful readjustment of properties of the time-even components of the 
mean fields. 



C. Modified gauge- violating Skyrme functionals 

In this Section we present results of calculations for Skyrme functionals with coupling 
constants which do not obey the gauge- invariance conditions (|2.10a]) - (|2.10cj) . In order to 



simplify the discussion, the gauge-invariant time-odd terms sf and s t -As t are disregarded, 
i.e., Ct=C^ s =0. Similarly as in the previous Section, we only consider simultaneous modi- 
fications of the isoscalar and isovector coupling constants. This leaves us with a possibility 
of modifying the time-odd coupling constants, C{, C t T , and while leaving unchanged 
their time-even gauge-invariance partners, C t r , C t J , and C^ J . 

First we have checked the separate role of the time-even and time-odd parts in the term 
SfT t —J\ which has been discussed in the previous Section. It turned out that the time-even 
part, J|, influences the results for rotational properties in a negligible way. Therefore, for 
all practical purposes the results discussed previously should be attributed to the time-odd 
part, s t -T t . This allows us to concentrate here on the two other time-odd terms, ji\ and 
Sf(Vxi t ). 

It is worth recalling at this point that in the gauge-invariant energy functional ([2.1 1|) 



the jl term is associated with the effective-mass term p t r t , and the term s t - (Vxj t ) comes 
together with the spin-orbit term p t (V-J t ). Therefore, the first one is absent in any theory 
which has the effective mass m* equal to the free mass m. This is the case for all studies 
based on phenomenological single-particle potentials, like the Nilsson or WS ones. 

On the other hand, the spin-orbit term is crucially important for a correct ordering of 
single-particle shells and is always taken into account in realistic calculations. However, the 
studies based on phenomenological potentials always disregard its time-odd gauge partner, 
and therefore should be considered as gauge-violating approaches, similar to the selfconsis- 
tent calculations considered in this Section. 

In Figs. H and |^ we present the results obtained for the same pairs of superdeformed 
bands as considered above, namely, for 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 151 Tb(l), respec- 
tively. Open symbols refer to the functionals with C{—0 and squares to those with C t Vjf =0. 
Consequently, the full circles repeat here the same results as those shown by full circles in 
Figs. || and f|. On the other hand, open squares correspond to the Skyrme functionals with 
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all time-odd terms neglected. 

As expected from the semiclassical analysis [43|, omitting the j 2 term leads to much 



smaller values of . The obtained decrease is of the order of 10-13%, depending on the 
value of the angular frequency, i.e., it is much smaller than the decrease of about 40% 
obtained from non-selfconsistent estimates [[EJ . With a very high precision, this decrease is 
identical for both nuclei in pairs differing by a hole in the 7r[301]l/2 (r=+i) orbital. Indeed, 
even if the values of change by as much as 13 /i 2 MeV _1 , the relative dynamical moments, 
Figs. ||(b) and ||(b), do not change at all, and are always well below 1 fr 2 MeV _1 . In addition, 
the j\ term has only a minor influence on the values of relative alignments 51. 

The presence of the second term studied here, Sf(Vxjt), has very little effect on the 
values and relative values of . However, its impact on the values of relative alignments 
is very large. Removing this term from the functional, decreases the values of 51 by about 
0.5 fr. For the 151 Tb(2)- 152 Dy(l) pair of bands, this leads to relative alignment which has 
the experimental value at the high-spin end of the band, Fig. |5](c). At the low-spin end 
the value of 51 is now underestimated by about 0.3 h. According to Eq. ( |3.1|) , this gives 
the differences of the 7-ray energies increasing from about — 6keV to 0, i.e., the identity 
of E 1 is still not reproduced with a sufficient precision. The sensitivity of results to the 
C t J ' coupling constant suggests, however, that a fit of this kind of term to experimental 
data may, in principle, improve the agreement. As mentioned previously, such a fit should 
take into account many different bands and cannot be based on a rather restricted set of 
examples discussed here. 

For the 150 Gd(2)- 151 Tb(l) pair of bands, the removal of the s t -(Vxj t ) term, while 
keeping the term j 2 unchanged, leads to a very good description of the values of J^ 2 \ and 
of the relative alignments 51 simultaneously, Fig. |6|. In this case, the relative alignments 
agree with the experiment to better than 0.1 7t, which corresponds to the identity of the 
7-ray energies to better than 1 keV. 



IV. CONCLUSIONS 

In the present study we have applied the Hartree-Fock cranking method with the Skyrme 
interactions to describe rotational states in selected superdeformed nuclei. We have analyzed 
in detail properties of two pairs of bands, namely, 151 Tb(2)- 152 Dy(l) and 150 Gd(2)- 151 Tb(l). 
Experimentally, these bands are pair- wise identical, i.e., the corresponding 7-transition en- 
ergies are identical within 2 keV each. 

In agreement with the previous interpretations, we have fixed the configurations of the 
151 Tb(l) band as a hole-structure created in the 7r[651]3/2 (r—+i) or 7164 Nilsson intruder 
orbital with respect to the magic superdeformed 1 |gDy 8 6 core. Similarly, the 151 Tb(2) and 
150 Gd(2) bands have been constructed as the 7r[301]l/2 (r=+i) hole configurations in the 
corresponding 152 Dy(l) and 151 Tb(l) cores. 

The dynamical moments calculated for the 151 Tb(2)- 152 Dy(l) pair of bands and for 
the 150 Gd(2)- 151 Tb(l) pair of bands are identical within 2h 2 MeV~ 1 . This result correctly 
reproduces the identity of these two pairs of bands obtained in experiment. It confirms 
that the internal structure of the 7r[301]l/2 (r=+i) orbital is responsible for the occurrence 
of these particular identical bands. The obtained identity of does not depend on the 
version of the Skyrme force used, neither it depends on including or disregarding various 
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time-odd components in the mean field of the rotating nucleus. The fact that the two bands 
calculated in 151 Tb, i.e., 151 Tb(l) and 151 Tb(2), are not identical (similarly as in experiment), 
indicates a crucial role of the single-particle structure of the involved orbitals, and shows 
that in the HF cranking theory the phenomenon of identical bands is not a generic, built- 
in-by- assumption result. 

The relative alignments of the studied pairs of bands are more difficult to reproduce than 
the simple identity of the dynamical moments. First of all, the results do depend on the 
version of Skyrme interaction and on the time-odd components included in the mean field. 
When the complete Skyrme functionals are used, i.e., when all time-odd terms are taken 
into account, the calculated relative alignments do not reproduce the experimental data. 
The disagreement obtained for the SkM* interaction is particularly large, while the SIII and 
SkP interactions also give too large relative alignments. Within the SkM* parametrization, 
good agreement with data is obtained when all time-odd terms are disregarded except of 
the one which involves the density of current, j, and is related to the effective-mass term by 
the gauge transformation. 

This specific result does not imply that the time-odd terms present in the formalism 
should be generally eliminated from its applications. By studying cases with some time-odd 
terms removed we only aimed at gaining information which may be useful when improving 
the parametrisation of effective interactions. 

As seen from the results obtained in the present study, a readjustment of parameters of 
effective forces in the time-odd channel seems to be necessary for a detailed description of 
the identical 7-ray transitions in superdeformed nuclei. However, this should be done by 
taking into account many available high-spin data simultaneously, and also should involve 
a readjustment of the time-even part (the spin-orbit interaction, in particular) which is 
responsible for the ordering of the single-particle energies and routhians. The work along 
these lines is now in progress. 
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APPENDIX A: COUPLING CONSTANTS OF THE SKYRME FUNCTIONAL 



In its standard form, the Skyrme interaction (see e.g. Ref. [12]) depends on 10 pa- 
rameters, t , x , t\, xi, t 2 , x 2 , £3, x 3 , a, and W. For xi=x 2 =0, the energy density which 
corresponds to the Skyrme interaction has been derived in Ref. 



26 . In terms of the isovec- 



tor and isoscalar coupling constants, and for the complete interaction, this energy density 
is given in Eqs. (p.2|)- 



The zero-order coupling constants, Cf and Cf, correspond to the velocity- independent 
terms of the interaction, and can be expressed by the to 5 ^ch t^, an d £3 parameters, as it 
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shows Table Il|. The zero-order time-odd coupling constants Cf (for t=0 and 1) are linear 
combinations of the zero-order time-even coupling constants Cf (Table |V|) . 

Similarly, the second-order coupling constants, C^ p , CJ, C t As , and Cf, correspond to 
the velocity-dependent terms of the interaction given by the parameters t X) X\, t 2 , and x 2 , 
and are presented in Table |V[ The second-order time-odd coupling constants, C^ s and 
C-f, are linear combinations of the second-order time-even coupling constants, C^ p and C\ 
(Table |VT|). Another four second-order coupling constants, Cf and Cl , which also depend on 
the same ti, x±, t 2 , and x 2 parameters, are given by the gauge- invariance conditions ([2.1 0a|) 
and P^UBQ . 

Finally, the second-order coupling constants C^ J are given by the spin-orbit term of the 
Skyrme interaction and depend on the parameter W, C^ J =— \W and CY J =—\W, while 
the other two, C t Vj , follow from the gauge- invariance condition ( ^.10c[ ). An extension of the 
Skyrme energy density, which introduces the coupling constants C^ J and CY J independent 
of one another, is discussed in Ref. fP 
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TABLES 



TABLE I. Coupling constants in the gauge-invariant energy density ( 2.11| ) obtained for the 
SIII, SkP, and SkM* Skyrme interactions. Zero-order coupling constants Cf and Cf ar given in 
MeV-fm 3 and the remaining (second-order) ones are given in MeV-fm 5 . 





t=0 


t=l 


SIII SkM* SkP 


SIII SkM* SkP 


Cf(p=0) 

C?(p=p nm ) 

ci 
cf J 


-423.3 -991.9 -1099.4 
-296.4 -237.7 -335.6 

-63.0 -68.2 -60.1 
44.4 34.7 0.0 

-90.0 -97.5 -75.0 


268.8 390.1 580.6 
141.2 150.8 188.4 
17.0 17.1 35.1 
-30.6 -34.1 -44.6 
-30.0 -32.5 -25.0 


C!( P =o) 

C t S 0=Pnm) 

c? s 
ci 


14.1 271.1 152.3 
56.4 31.7 -31.4 
17.0 17.1 -4.2 
-30.6 -34.1 7.7 


141.1 330.6 366.5 
98.8 91.2 78.5 
17.0 17.1 9.8 

-30.6 -34.1 -41.1 



TABLE II. Proton quadrupole moments (in barns) calculated for the SIII, SkM*, and SkP 



Skyrme forces at uj=0.5M.eV /h. 





SIII 


SkM* 


SkP 


152 Dy(l) 


18.55 


18.25 


18.36 


151 Tb(2) 


18.69 


18.38 


18.48 


151 Tb(l) 


17.59 


17.22 


17.41 


150 Gd(2) 


17.74 


17.36 


17.54 



TABLE III. Zero-order coupling constants as functions of parameters of the Skyrme interac- 
tions, expressed by the formula: C=^(ato+btoXo)+-^p a (at3+bt^X3). 





a 


b 


c p . 


3 





ci 


-1 


-2 


Co 


-1 


2 


Cf 


-1 






TABLE IV. Zero-order time-odd coupling constants as functions of the time-even coupling 
constants, expressed by the formula: C=|(aCQ+6Cf ). 





a 


b 




-2 


-3 


Cf 


-1 
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TABLE V. Second-order coupling constants as functions of parameters of the Skyrme interac- 
tions, expressed by the formula: C=^(ati+btiXi+ct2+dt2X2)- 





a 


b 


c d 


Co 

c\ 


-9 
3 
12 

-4 



6 


-8 


5 4 
1 2 
20 16 

4 8 


r<As 
°o 


3 
3 
-4 
-4 


-6 

8 



1 2 
1 
4 8 
4 


TABLE VI. Second-order time-odd coupling constants as functions of the time-even coupling 
constants, expressed by the formula: C=^(aC^ p +bC^ p +cC$ +dC{) . 




a 


b 


c d 


cl 
cl 



-4 


16 


6 
-4 
48 
■16 


3 9 

3 -3 
-4 12 

4 -12 
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FIGURES 



FIG. 1. Calculated dynamical moment for the yrast band of 152 Dy, part (a), the relative 
dynamical moment calculated for the 151 Tb(2) and 152 Dy(l) bands, part (b), and the relative 

alignment 51 between these two bands, part (c). The experimental points are denoted by asterisks. 
Complete Skyrme functionals of SIII, SkM*, and SkP interactions have been used. Note the scale 
in (b) expanded five times as compared to (a). 

FIG. 2. Same as Fig. [l], but for the 151 Tb(l) band, (a), and for the differences between the 
150 Gd(2) and 151 Tb(l) bands, (b) and (c). 

FIG. 3. Same as Fig. [I], but for the complete SkM* functional (open squares) compared to 
three versions of the modified SkM* functional, namely, C^=0 (open circles), Cf=C^ s =0 (full 
squares), and Cj =C^ s =Cf=Q (full circles). The modified functionals are gauge invariant, i.e., 
C/=0 whenever Cf =0, Eq. (pT0b| ). 

FIG. 4. Same as Fig. |, but for the 150 Gd(2)- 151 Tb(l) pair of bands. 

FIG. 5. Same as Fig. ||, but for the modified gauge-invariant SkM* functional given by 
C^=C^ s =C*=0 (full circles), compared to three versions of the gauge- violating SkM* functional 
obtained by putting in addition C{=0 (open symbols) and/or C t 3 =0 (squares). The version with 
C/=C t V,7 =0 (open squares) corresponds to all mean- field time-odd terms neglected. 

FIG. 6. Same as Fig. g but for the 150 Gd(2)- 151 Tb(l) pair of bands. 
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Cranked Hartree-Fock - Skyrme SkM* 
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